Realization of an effective ultrahigh magnetic field on a nanoscale 
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In tunnel junctions of which at least one side is a ferromagnet, very large magnetic polarization 
change (~ O.lpn) and splitting of the spin up and spin down Fermi energy (~ O.leV) can be created 
under steady state finite current conditions (bias voltage ~ 1 volt). This is much higher than can 
be created by the highest magnetic field on earth. We illustrate this with a specific calculation of 
a recently observed very large Hall effect in the Al side of a Co-I-Al tunnel junction. Other recent 
experiments that support this idea are discussed. 

I. INTRODUCTION AND SUMMARY 

There is recently much interest in spin polarized transport in ferromagnetic tunnel junctions, motivated by the high 
D.C. magnetoresistance observed, Jl],p| and their potential applications as magnetic sensors and nonvolatile memories. 
These structures are junctions of two metals with a thin (« 10A ) layer of an insulator in between. They are 
symbolically denoted as an F-I-F structure where I stands for insulator. The spin polarized transport in this system 
has been studied theoretically without taking into account the electron-electron (e-e) interation [|-|| and including 
the e-e interaction. @-|| It is found that in tunnel junctions at least one side of which is a ferromagnet, under 
steady state moderate currents (« 10 -5 Ampere) there will be a very large change in the polarization (ss O.lps) 
and splitting of the spin up and spin down Fermi energies (« O.leV). 

The implication of this splitting has not been explored. By applying different voltages, the Fermi levels can be 
moved and different parts of the band structure can be probed. Different physical characteristics of different parts of 
the band structure can be manifested. One can envision inducing different phase transitions by tuning the applied 
voltage. 

Experimental results that illustrate the new physics that can happen as a result of the splitting are beginning 
to appear. The big splitting of the spin up and spin down Fermi energies provides for an explanation of the bias 
dependence that is experimentally observed. || We discussed here an illustration of the probing of a different part 
of the band structure with a detailed calculation of a giant anomolous Hall resistivity recently discovered in Co-I-Al 
tunnel junctions. jl0| Good agreements with the experimental results are obtained. We then conclude by discussing 
two other experiments that provide further illustration for the conclusion discussed here. 

The changes in the Fermi levels can in principle also be produced by external magnetic fields of strength of the 
order of 1000 Teslas. This field strength is much higher than can be created by the highest magnetic field on earth. 
This presents a possibility of a new area of research of ultrahigh magnetic field physics. 

That a current can induce a magnetization in a ferromagnetic- paramagnetic junction (no insulator in between) 
was first discussed by Johnson and Silsbee || and van Son et al. j|. The magnetized current from the ferromagnet 
induces a magnetization in the paramagnet due to a spin bottleneck effect and causes a splitting of the spin up and 
the spin down Fermi energies A/i. The ratio of this splitting to the driving current / is of the order of the resistance 
of the metal. This same reasoning can also be applied to the F-I-P structure with an insulator in between, and 
the magnetization and the splitting is of the same order of magnitude. The resulting magnetization from this "spin 
accumulation effect" is very small and is different from the physics discussed here. For the F-I-P structure and our 
mechanism, the ratio Ap/I is of the order of the resistance of the insulating barrier. This is several orders of 
magnitude larger than that expected from the spin accumulation effect. 

II. UNDERLYING FRAMEWORK 

We first recapitulate briefly the reasoning that led to the splitting of the Fermi levels. The physics can be approxi- 
matly described by the following sets of equations. The first is that of global charge current conservation: 

V • J = -dp/dt (1) 

Here J, p are the total (including both spins) current and charge densities. The second equation is the diffusion 
equation which expresses the current density for spin s, J s , as a sum of (a) the external driving current Jq s = a s E 
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that is controlled by the conductivity a s and the external electric field E; (b) the current driven by a density gradient 
(diffusion) expressed as a gradient of the chemical potential V// s = Vp s /N s through the density of states at the Fermi 
surface N s and (c) the current driven by the internal electric "screening" field VWo: |T^ ] 

J s =ct s [V/x s - VW + eE]/e (2) 

Here a s is the conductivity for spin channel s. Wq(t) is the local electric (screening) potential due to the other electric 
charges that is determined self-consistently: Wo(r) = J d 3 r'U(r — r')p(r') where U is the Coulomb potential. Because 
of this self-consistent screening, charge fluctuation dies off with a length scale A called the screening length, which 
is of the order of 10A for metals. Finally, the magnetization density changes M relaxes with a length scale of l s f 
where the bar indicates a renormalization effect due to the electric field: 

V 2 M - M/P sf = (3) 

For the systems of interest, the spin diffusion length l s f, is of the order of 1000 to 10000 A and is much larger than 
the screening length A. Because of these two very different length scales, there is a delicate balance in the current 
from the charge and the magnetization diffusion that has to be maintained. 

One can solve these sets of equations on each side of the junction. The currents on opposite sides of the junction 
are related by the boundary condition that the current J s is continuous across the junction and that 

Afi s - AW = r(l - s 7 ) J s (4) 

where r(l — sj) is the resistivity of the junction for spin s = ±1. 

The final solution can be developed in increasing powers of X/l s f- Solving equations (1-3) for a junction of thickness 
d at the origin, we find "dipole layers" of charge (p) and magnetization (M) densities that peaks at the junction and 
dies off exponentially away from the junction with a functional form /(z, /) = cxp[— (|z| — d/2)/l]: p"3| ] 

p « Pl f(z, \)X/l sf + p 2 f(z, l s f)X 2 /P sf 



M « M f(z, l s f). 

In this paper, we use units of M in terms of the Bohr magneton. The charge and magnetization densities are now 
coupled: pw = Gp2o where G is a coefficient of the order of unity. |l3 p 2 o — MqDm/Dd- Here Dm {Dd) is the 



diffusion constant for the magnetization (charge) that can be expressed in terms of the conductivities. |16| The charge 
dipole layer densities are smaller than the magnetization densities by a factor \/l s f*. This comes from the two very 
different length scale of change for the charge and the magnetization degrees of freedom. There is a small correction 
to the charge density that is of much longer range than the screening length. In general the constants pi, Mq on 
opposite sides of the junction can be different from each other. The total charge on opposite sides of the junction 
are opposite in sign and equal in magnitude, of course. 

After matching the boundary condition (eq. (4)) one obtains for the magnetization and the Fermi energy splitting: 



Mo = FrJ (5) 



p+ - p-= Q.5M{1/N + + 1/AT_) (6) 

where F is a coefficient of the order of the inverse density of states at the Fermi energy |l8| . It is different for parallel 
and antiparallel alignment of the magnetizations on opposite sides of the junction. Thus the magnetization density 
induced is controlled by the resistance of the junction, not of the metal! To illustarate the application of this, we 
first calculate the Hall conductivity <j xy of Al in a Co-I-Al junction. 



III. ANOMLOUS GIANT HALL RESISTIVITY 

Otani et al. [0 recently measured the Hall coefficient of Al in a Co-I-Al structure. They found that above a 
threshold voltage of about 1 volt, there is a rapid increase in the Hall voltage that is of the same sign (negative) 
as but much larger than the ordinary Hall coefficient of Al. This effect disappears as the temperature is increased 
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from 30K to 150 K. Their original explanation in terms of an extraordinary Hall effect is inapplicable with a revised 
estimate of the Hall resistivity that is several orders of magnitude higher. 

In our picture, as the external voltage V is increased, the spin up and spin down Fermi energies will move apart by 
an amount 2 A proportional to V. The sample is polycrystalline experimentally (To), Without detail knowledge of the 
nature of this polycrystalline sample, for simplicity we consider the crystal structure as a single fee crystal of bulk Al 
with lattice parameter 4.049lA, and have taken the direction of the magnetization parallel to the z symmetry axis of 
Al. We have calculated the Hall conductivity a xy from first principles as a function of A. 

Our calculations are performed using the band structure obtained with the self-consistent full-potential linearized 
augmented-plane-wave (FLAPW) method jl9| under the generalized gradient approximation (GGA) [^0| with spin- 
orbit coupling [pl|-p3|. The conductivity is calculated using the Kubo formula. The Brillouin zone sampling is 
performed using 4000 special k-point meshes, which yielded 315 points in the irreducible Brillouin zone. j24| The 
dependence of conductivity on the frequency have been tested on the Kerr effect of bec Fe. For an inverse lifetime 
r = 0.45 ~ Q.QbeV , our result is in very good agreement with the experimental results. The Hall conductivity is the 
zero frequenct limit of the conductivity. 

Our result for the Hall conductivity is shown in Fig. 1. The sign and the shape (including the threshold) are in 
good agreement with experiment. The physics of the threshold and the temperature dependence can be understood 
as follows. The band structure of Al with the spin up and spin down Fermi levels is shown in Fig 2 for the case 
A = 0.2eV. Two band crosses along the direction XZW. The crossing is spin-orbit split. At zero magnetization, 
this corssing is below the Fermi energy. When the spin-up and the spin down Fermi levels move apart, one of them 
will move towards this crossing and a significant interband contribution will come in for the Hall conductivity. The 
thresold is determined by the distance between the energy at this crossing and the Fermi energy at zero field. Because 
the spin-orbit band gap is very small, we expect the interband contribution will decrease significantly when the 
temperature becomes comparable to this band gap. 

The magnitude of our result depends on the scattering time r = Ti/T. Our result is in the range of 1 to 10 1/Q — cm. 
Because of the nonuniformity of the current flow due to the geometry of the experimental sample, and becuase we 
do not know the detail material characteristic of the polycrystalline experimental sample, we do not expect to get 
exact quantitative agreement with the experimental results. If we pick a typical T = 0.25eV, our result is in good 
agreement with experimentally revised value of p xy — 1.5 x 10 _8 f2 — m for an average p xx = 3.4 x 10 _6 f2 m. (Recall 
that a xy = p xy I ' p 2 xx ) This diagonal resistivity is consistent with recent unpublished results by the same authors for 
a ohmic (F-P) Co-Al junction. 

Our picture can be further tested experimentally. In our calculation, the Hall resistance goes through a maximum 
and decreases as the bias voltage is further increased. So far the experiment was carried out up to a bias when the 
maximum of p xy was reached. We predict that the Hall resistivity will decrease as the voltage is further increased. 
We next discuss some other recent experimental results that supports the present picture. 

Moodera and coworkers recently found a fascinating series of bias dependences, including a negative junction 
magnetoresistance, for a F-I-P-F structure with varying thickness for the paramagnet. |2(| To explain their result, 
they assuumed a model with different spin up and spin down Fermi energies. The result discussed here provides for 
a possible justification of the model that they have used. 

Ono, Shimada and Ootuka |27|] recently studied the magnetoresistance of ferromagnetic single electron transistors 
made from tunnel structures. Depending on the gate voltage, this device can be in an on or an off state. In the on- 
state, they found a magnetoresistance ratio of about 4 %. They observed a ten fold increase of the magnetoresistance 
ratio in the off-state. In the off-state, the transport depends exponentially on the activation barrier which is a sum 
of the Coulomb charging energy and the difference of the Fermi energy on the left and the right hand side of the 
junction. ]28[ ] The large magnetoresistance can be understood in the present picture from the different splitting of the 
Fermi levels between the cases of the parallel and the antiparallel alignment of the magnetizations. 

In conclusion, we propose that magnetic tunnel junction provides for new opportunities to observe high magnetic 
field physics. Theoretical and experimental evidences are provided to support this claim. 
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FIG. 1. Dependence of conductivity and the spin polarization on the spin-up spin-down Fermi level splitting 2A = E F p — E P n . 

FIG. 2. Energy bands for fee Al in the ferromagnetic ground state with Fermi level splitting 2A = 0.4eV. The Fermi levels 
are shown by the solid horizontal lines. 
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